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RESEARCH  AMBIENT  NOISE  DIRECTIONALITY  (RANDI)  3.1 
PHYSICS  DESCRIPTION 


1.0  INTRODUCTION 

The  purpose  of  this  report  is  to  present  a  description  of  the  physics  underlying  the  Research 
Ambient  Noise  Directionality  noise  model  version  3.1,  hereafter  referred  to  as  RANDI  3.1.  A 
companion  report,  the  RANDI  3.1  Users  Guide  (Breeding  et  al.  1994),  has  been  prepared  for  those 
who  plan  to  use  this  model. 

The  RANDI  3.1  noise  model  is  designed  to  predict  the  response  of  low-  to  mid-frequency  sonar 
receivers  to  the  ocean  acoustic  noise  field  in  locations  with  highly  variable  surrounding  bathymetry 
and  range-dependent  sound  speed  structure.  Such  environments  are  typical  of  the  complex  oceanographic 
areas  found  in  shallow  water  and  coastal  areas.  However,  the  RANDI  3.1  model  can  be  applied 
equally  well  in  deep-water  areas.  In  RANDI  3.1,  state-of-the-art  acoustic  propagation  models, 
computational  algorithms,  and  geographic  environmental  data  bases  are  combined  into  a  single  product. 

The  RANDI  3.1  model  is  based  on  the  Research  Ambient  Noise  Directionality  II  (RANDI  II) 
model  developed  at  the  SACLANT  Undersea  Research  Centre  by  Hamson  and  Wagstaff  (1983). 
The  RANDI  II  model  is  an  ambient  noise  model  that  predicts  noise  levels  and  directionalities  for 
user-specified  environmental  and  shipping  conditions  using  adiabatic  mode  theory  to  propagate 
energy  from  individual  ships  to  a  receiver  array.  The  array  response  for  each  ship  is  calculated  by 
summing  the  complex  pressure  due  to  individual  modes  at  each  hydrophone  in  the  receiver  array, 
followed  by  either  coherent  or  incoherent  summation  across  modes. 

The  RANDI  II  model  suffered  from  computational  limitations  centered  around  the  use  of  normal 
mode  theory,  making  it  difficult  to  use  the  model  in  deep  water  where  the  number  of  propagating  modes 
is  large  or  in  areas  with  rapid  oceanographic  or  bathymetric  changes,  such  as  found  in  shallow 
water.  The  RANDI  3.1  model  overcomes  these  difficulties  by  using  parabolic  equation  (PE)  propagation 
loss  models  to  propagate  energy  from  individual  ships  to  the  receiver  array.  Currently,  RANDI  3.1 
incorporates  two  PE  models:  (1)  finite  element  PE  (FEPE)  and  (2)  the  Navy-Standard  PE  (NSPE), 
also  referred  to  as  the  split-step  PE  (SSPE).  The  SSPE  model  (Hardin  and  Tapper!  1973)  is  a 
solution  to  the  parabolic  equation  that  involves  a  marching  type  Fourier  transform.  The  SSPE  has 
a  relatively  wide  angle  capability  with  up  to  40°  half-beamwidth.  The  FEPE  model  (Collins  1988, 
1989)  is  based  on  a  Pade  series  expansion  of  the  depth-dependent  part  of  the  parabolic  equation. 
The  FEPE  was  developed  in  an  effort  to  handle  propagation  angles  approaching  90°.  The  use  of  PE 
propagation  codes  in  RANDI  3.1  greatly  extends  the  frequency  range  and  the  range  of  oceanographic 
conditions  for  which  valid  noise  predictions  can  be  made. 

The  RANDI  3.1  model  computes  the  shipping  noise  complex  pressure  by  one  of  three  propagation 
methods.  These  are  the  rigorous  method,  the  radial  accumulation  method,  and  the  high-resolution 
radial  accumulation  method.  The  rigorous  method  executes  the  propagation  model,  either  FEPE  or 
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SSPE,  from  each  ship  to  the  receiver.  The  radial  accumulation  method  discretizes  the  azimuthal 
direction  into  radials  and  executes  the  propagation  model  from  the  maximum  ship  range  of  each 
radial  to  the  receiver,  accumulating  ships  as  it  marches  inward.  The  high-resolution  radial  accumulation 
method,  in  addition  to  executing  the  radial  accumulation  method,  computes  the  complex  pressures 
at  each  hydrophone  for  each  ship  instead  of  for  each  radial.  The  radial  accumulation  method  and 
the  high-resolution  radial  method  are  computed  using  only  the  FEPE  model. 

In  RANDI  3.1,  as  well  as  RANDI II,  local  wind  noise  is  estimated  using  a  wave  theory  model 
developed  by  Kuperman  and  Ingenito  (1980).  A  cross  spectral  density  matrix  is  determined  for  all 
hydrophones  in  the  receiver  array. 

The  RANDI  3.1  model  is  fully  interfaced  to  supporting  geographical  environmental  data  bases 
of  bathymetry,  sound  speed,  bottom  properties,  shipping,  and  wind  speed.  The  model  extracts 
environmental  information  from  Navy-Standard  data  bases  including  Earth  TOPOgraphy  5  (ET0P05) 
and  Digital  Bathymetry  Data  Base  Confidential  (DBDBC)  for  bathymetry;  the  Historical  Ocean 
Profiles  (HOP),  the  Provinced  Generalized  Digital  Environmental  Model  (Provinced  GDEM),  and 
the  High  Resolution  Shallow  Water  Sound  Speed  (SWSS)  data  bases  for  sound  speed  profiles; 
and  the  Low  Frequency  Bottom  Loss  (LFBL)  and  Consolidated  Bottom  Loss  Upgrade  (BLUG)  data 
bases  for  bottom  parameters.  Shipping  information  is  obtained  from  the  Historical  Interim  Temporal 
Shipping  (HITS)  3.0  and  3.1  data  bases,  and  wind  speed  from  the  Historical  Wind  Speed  (HWS) 
data  base.  The  RANDI  3.1  model  also  allows  for  a  variety  of  array  types,  beamforming  options,  and 
output  displays. 

The  shipping  noise  sources  are  described  in  Sec.  2.0  and  the  environmental  data  bases  are  described 
in  Sec.  3.0.  In  Sec.  4.0,  the  shipping  noise  propagation  models  are  explained  and  the  wind  noise 
model  is  presented  in  Sec.  5.0.  Noise  predictions  that  are  independent  of  the  receiver  are  discussed 
in  Sec.  6.0,  while  noise  predictions  that  depend  upon  the  receiver  array  are  found  in  Sec.  7.0. 


2.0  SHIPPING  NOISE  SOURCES 

By  far,  the  major  component  of  the  total  low-frequency  (<500  Hz)  ambient  noise  field  is 
the  noise  generated  by  ship  traffic.  Even  very  distant  shipping  can  be  the  dominant  noise  near  the 
horizontal  at  a  receiver  array  if  there  is  a  shoaling  sound  channel  axis  or  downslope  conversion 
(Wagstaff  1981).  Before  it  is  possible  to  determine  the  noise  due  to  shipping  at  a  receiver,  it  is 
necessary  to  determine  where  the  ships  are  located.  Each  ship  is  then  assigned  a  heading  and  source 
level,  as  is  explained  below. 


2.1  Estimation  of  Ship  Locations  and  Headings  from  HITS 

The  most  extensive  data  bases  of  shipping  available  are  the  Historical  Temporal  Shipping 
(HITS  3.0  and  3.1)  data  bases  (Molinelli  1990;  Naval  Oceanographic  Office  1993;  Barnes  1995), 
which  contain  shipping  densities.  These  data  bases  have  been  incorporated  into  the  RANDI  3.1 
model.  However,  the  RANDI  3.1  model  requires  discrete  ship  locations  to  calculate  shipping  noise, 
so  the  densities  are  converted  to  discrete  locations. 

The  HITS  data  bases  provide  the  shipping  density  for  each  1®  grid.  The  densities  contain  the 
expected  number  of  ships  within  each  grid  for  a  month,  season,  or  year.  Five  types  of  ships  are 
identified  according  to  size:  merchant  ships,  tankers,  large  tankers,  super  tankers,  and  fishing 
vessels.  The  fishing  vessel  densities  are  given  in  the  HITS  3.0  data  base,  and  they  are  averaged 
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annually  and  for  the  winter  and  summer  seasons.  The  densities  for  the  merchant  ships,  tankers, 
large  tankers,  and  super  tankers  are  listed  in  the  HITS  3.1  data  base,  and  they  are  averaged  seasonally, 
monthly,  and  yearly.  Both  data  bases  are  unclassified.  Figure  1  gives  an  example  of  the  shipping 
densities  extracted  within  approximately  3000  km  of  a  receiver  array  located  in  the  East  China  Sea 
outside  the  Korea  Strait  at  31°  N,  127°  E.  The  densities  shown  are  the  sum  of  densities  for  the  five 
individual  ship  classes  included  in  the  HITS  data  bases.  Densities  for  each  ship  class  can  also  be 
extracted.  For  example,  Fig.  2  depicts  the  shipping  densities  for  fishing  vessels  in  this  area. 

For  a  given  ship  type,  the  shipping  densities  for  a  region  are  used  to  compute  the  discrete 
number  of  ships  for  each  1°  grid.  This  is  done  using  a  Poisson  random  distribution  function  with 
the  extracted  ship  density  as  the  mean.  The  Poisson  equation  is 


ships 

^  e  ^\x 
r<  Z 


k=\ 


k\ 


(1) 


where  r  is  a  uniformly  distributed  random  number  between  0.025  and  0.975,  p.  is  the  shipping 
density,  and  ships  is  the  number  of  discrete  ships.  The  above  equation  is  solved  with  the  condition 
ships  =  0  if  r  <  e~^.  The  random  number  r  is  confined  to  the  range  [0.025,  0.975]  rather  than  [0,1] 
to  prevent  the  noise  field  from  being  dominated  by  extremes  of  the  Poisson  distribution. 


Fig,  1  —  HITS  3.1  shipping  densities  in  the  East  China  Sea.  All  five  classes  of  vessels  are  included. 
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Fig.  2  —  HITS  3.1  fishing  vessel  densities  in  the  East  China  Sea 


For  each  discrete  ship  found  in  the  grid,  an  exact  ship  location  is  computed  within  that  grid  cell 
using  a  uniform  distribution  function.  If  a  grid  cell  contains  both  land  and  water,  the  ship  location 
is  adjusted  so  as  to  not  lie  on  land.  Also,  the  great  circle  path  distance  and  bearing  from  the  ship 
location  to  the  receiver  is  calculated.  Figure  3  displays  the  discrete  shipping  positions  found  for  the 
shipping  densities  shown  in  Fig.  1.  Each  ship  type  is  represented  by  a  unique  symbol.  Ships  or 
radials  that  are  blocked  from  the  receiver  by  land  are  ignored  in  the  calculation  of  noise. 

With  the  exception  of  fishing  vessels,  a  ship  heading  is  computed  for  each  ship  by  determining 
a  direction  perpendicular  to  the  gradient  of  the  shipping  density  in  the  region  defined  by  the  grid 
containing  the  ship  and  the  grids  immediately  surrounding  this  grid.  Assuming  the  gradient  points 
towards  the  shipping  lane,  the  direction  perpendicular  to  the  gradient  causes  the  ship  to  move 
parallel  to  the  shipping  lane  if  the  ship  is  dead-reckoned.  The  fishing  vessels  are  assigned  headings 
generated  from  a  uniform  distribution  function  to  simulate  their  meandering  nature. 


2.2  Ship  Source  Levels 

The  source  level  of  each  ship  depends  upon  the  length  and  speed  of  the  ship.  A  ship  length  and 
speed  are  calculated  for  each  ship  by  using  a  uniform  random  distribution  function,  but  forcing  the 
length  and  speed  to  have  the  limits  shown  in  Table  1. 
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Fig.  3  —  Discrete  ships  in  the  East  China  Sea  calculated  with  a  Poisson  distribution 


Table  1  —  Length  and  Speed  Limits 


SHIP  TYPE 

LENGTH  (ft) 

SPEED  (kt) 

Fishing  Vessel 

50-150 

7-10 

Merchant 

275-400 

10-15 

Tanker 

400-500 

12-16 

Large  Tanker 

500-700 

15-18 

Super  Tanker 

800-1200 

15-22 

The  ship  source  levels  are  computed  by  defining  an  “average”  ship  as  one  with  a  speed  of  12  kt 
and  a  length  of  300  ft.  This  average  ship  is  assigned  a  source  level  of  L^o  dB  as  a  function  of 
frequency.  The  noise  source  levels  of  the  actual  ships  in  the  model  are  then  calculated,  on  the  basis 
of  their  individual  speeds  and  lengths,  by  the  following  empirical  equation  based  on  Ross  (1987). 

Ls  (f,  V,  Is)  =  Lso  (/)  +  60  log(v/12)  +  20  logC^/SOO)  +  df*dl  + 3.0,  (2) 

where  v  is  the  ship’s  speed  in  knots  and  4  is  its  length  in  feet  and 
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I  0.00  </<  28.40  #=8.1 

^  I  28.4  </<  191.6  df  =22.3-9.11  *  log  (f)  ^ 

dl  =  3643.0. 

For  convenience,  the  equation  for  Lso  if)  when  /<  500  Hz  is  given  by 

i^io(/)  =  -10  log(10-106  log/-14.34^  J03.32  log/ -21.425)  (4) 

For  values  of  frequency  greater  than  500  Hz,  Lso  if)  is  given  by 

4,  (/)=  173.2 -18.0  log(/).  (5) 

Table  2  shows  representative  source  levels  for  the  different  types  of  ships  for  several  frequencies 
using  average  lengths  and  speeds  from  Table  1.  The  source  levels  are  in  units  of  decibels  referenced 
to  1  [0,  Pa. 


3.0  ENVIRONMENTAL  DATA 
BASES 

To  determine  the  amount  of  energy 
that  propagates  from  sources  to  receiv¬ 
ers,  it  is  necessary  to  characterize  the 
different  environments  with  which 
the  waves  (rays)  interact.  The  RANDI  3.1 
noise  model  is  designed  to  automatically 
extract  environmental  data  from  a  number 
of  data  bases  along  a  great  circle  path 
from  a  receiver  location  to  a  source.  A 
great  circle  path  is  the  shortest  distance  between  the  receiver  and  the  source  on  a  spherical  surface. 
The  great  circle  path  distance  {GCD)  from  the  receiver  to  the  source  is  given  by 

GCD  =  60  cos“^  (sin05  sin0^  +  cos0j  cosQr  cos(<I>^  -  4)^)) ,  (6) 

where  GCD  is  in  nautical  miles,  0^  is  the  source  latitude,  0;-  is  the  receiver  latitude,  is  the 
source  longitude,  and  O;-  is  the  receiver  longitude. 

Along  each  great  circle  path,  extraction  routines  are  used  to  retrieve  bathymetry,  sound  speed, 
and  bottom  characteristics  at  each  range  segment.  The  environmental  data  are  extracted  whenever 
the  great  circle  path  enters  a  new  grid  cell  of  a  particular  data  base.  This  section  gives  a  brief 
description  of  the  individual  data  bases  which  RANDI  3.1  accesses  to  simulate  the  propagation 
environment  as  accurately  as  possible. 


Table  2  —  Representative  Source  Levels 


SHIP  TYPE 

10  Hz 

25  Hz 

50  Hz 

100  Hz 

300  Hz 

Fishing  Vessel 

142.7 

146.5 

144.8 

136.0 

120.0 

Merchant  Ship 

160.9 

167.8 

162.6 

153.5 

137.1 

Tanker 

167.0 

170.8 

168.6 

159.2 

141.6 

Large  Tanker 

174.8 

178.6 

176.0 

166.3 

149.3 

Super  Tanker 

185.0 

188.8 

185.4 

174.6 

156.8 

3.1  Bathymetry 

The  Digital  Bathymetry  Data  Base  (DBDB)  (Naval  Oceanographic  Office  1987)  contains 
acoustically  measured  ocean  depths,  referenced  to  a  1500  m/s  sound  speed.  Two  DBDB  data  bases 


RANDI  3. 1  Physics  Description 


7 


are  incorporated  in  the  RANDI  3.1  model.  One  version  is  an  unclassified  bathymetry  data  base 
referred  to  as  ET0P05,  and  the  other  version  is  a  confidential  bathymetry  data  base  referred  to  as 
DBDBC.  Both  data  bases  contain  worldwide  bathymetry  data  measured  in  meters  at  a  1/12°  resolution. 
Figure  4  shows  a  view  of  the  ET0P05  bathymetry  found  in  the  East  China  Sea  outside  the  Korea  Strait. 


3.2  Sound  Speed 

The  RANDI  3.1  model  provides  two  data  bases  for  the  great  circle  path  extractions  of  sound 
speed.  One  of  the  data  bases  is  the  Provinced  GDEM  data  base,  which  is  also  known  as  the  HOP 
data  base  (Naval  Oceanographic  Office  1990b).  This  is  a  provinced  subset  of  the  GDEM  data  base. 
The  HOP  data  base  contains  temperature,  salinity,  and  sound  speed  profiles  at  standard  depths  in 
the  Atlantic  Ocean,  Pacific  Ocean,  Indian  Ocean,  and  Mediterranean  Sea.  The  standard  depths  are 
0,  10,  20,  30,  50,  75,  100,  125,  150,  200,  250,  300,  400,  500,  600,  700,  800,  900,  1000,  1100,  1200, 
1300,  1400,  1500,  1750,  2000,  2500,  3000,  4000,  5000,  6000,  7000,  8000,  and  9000  m.  The  HOP 
data  base  is  unclassified,  and  the  data  base  profiles  are  seasonal  with  a  resolution  of  1/2°. 

The  second  data  base  provided  in  RANDI  3.1  is  the  SWSS  data  base,  and  it  describes  the 
seasonal  temperatures  and  salinity  profiles  in  shallow  water  at  the  standard  depths.  This  data  base 
is  detailed  in  Crout  (1991)  and  was  developed  by  both  the  Office  of  Naval  Research  and  Air 
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Fig.  4  —  Bathymetry  in  the  East  China  Sea  as  extracted  from  the  ETOP05  data  base 
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Defense  Initiative.  The  SWSS  data  base  was  developed  using  the  GDEM  and  Master  Oceanographic 
Observation  Data  Set  data  bases  in  computing  the  shallow-water  profiles  for  depths  between  50  to 
500  m.  The  purpose  of  the  SWSS  data  base  is  to  allow  a  smooth  transition  between  the  GDEM  data 
base  profiles  and  shallow-water  profiles.  The  SWSS  data  base  has  a  1/12°  resolution,  is  seasonal, 
and  unclassified.  The  coverage  is  Northern  Hemisphere  shallow-water  regions,  excluding  the  Arctic 
region. 

An  example  of  sound  profiles  extracted  for  a  tract  near  the  Pacific  Coast  of  the  United  States 
is  shown  in  Fig.  5.  The  bathymetry  and  the  sound  channel  axis  are  also  shown. 

3.3  Bottom  Characteristics 

The  bottom  characteristics  data  bases  incorporated  in  the  RANDI  3.1  model  include  the  LFBL 
and  Consolidated  BLUG  data  bases.  The  LFBL  data  base  is  described  by  the  Naval  Oceanographic 
Office  (1990)  and  by  Spofford  et  al.  (1983).  This  data  base  contains  the  acoustic  properties  of  the 
ocean  sediment  and  the  sediment  thickness  at  low  frequencies  (50  to  1600  Hz).  In  the  data  base, 
the  ocean  is  divided  into  geographical  areas  that  are  homogeneous  in  sediment  and  composition 
type.  These  provinces  point  to  the  geoacoustic  parameters  which  provide  the  BLUG  geoacoustic 


Fig.  5  —  Bathymetry  and  sound  speed  profiles  for  a  tract  near  the  Pacific  Coast  of  the  United  States 
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profile.  The  geoacoustic  parameters  describe  the  sound  interaction  with  the  seafloor,  density,  sound 
speed,  attenuation  of  the  fluid  layer,  surface  layer  thickness  and  density,  and  substrate  reflectivity. 
In  addition,  for  shallow-water  provinces,  an  attenuation  exponent  is  included  as  a  geoacoustic 
parameter.  The  LFBL  data  base  is  classified  and  covers  the  Northern  Hemisphere,  the  Indian  Ocean 
to  50°  S,  and  the  South  Atlantic  Ocean  to  50°  S  at  1/12°  intervals. 

The  Consolidated  BLUG  data  base  is  an  unclassified  version  of  the  LFBL  data  base.  The 
Consolidated  BLUG  data  base  contains  the  same  parameters  as  the  LFBL  data  base  with  a  resolution 
of  1/12°  grid  for  the  Northern  Hemisphere  to  65°  N,  the  Atlantic  Ocean  to  70°  N,  and  the  Indian 
Ocean  to  10°  S. 

The  geoacoustic  parameters  contained  in  both  the  LFBL  and  the  Consolidated  BLUG  data  bases 
are  listed  in  Table  3  and  are  displayed  in  Fig.  6. 

Table  3  —  Geoacoustic  Parameters 

Sediment  Thickness  (s) 

Water/Sediment  Sound  Speed  Ratio  (pure  number) 

Thin  Layer  Thickness  (m) 

Thin  Layer  Density  (gm/cc) 

Sediment  Surface  Density  (gm/cc) 

Initial  Sound  Speed  Gradient  in  Sediment  (1/s) 

Sediment  Sound  Speed  Profile  Curvature  (1/s) 

Attenuation  at  z  =  0  (dB/m/kHz) 

Attenuation  Gradient  (dB/m/kHz/m) 

Basement  Reflection  Coefficient  (pure  number) 

Attenuation  Exponent  for  Shallow  Water  (pure  number) 

Two-Way  Travel  Time  (s) 
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Fig.  6  —  Bottom  parameters  included  in  the  LFBL  and  consolidated  BLUG  data  bases 
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Once  the  bottom  characteristics  are  extracted  from  a  data  base  they  are  then  converted  into  a 
sediment  sound  speed  profile  and  are  used  to  compute  the  density  and  attenuation  in  the  sediment 
and  bottom.  A  shear  sound  speed  and  attenuation  are  also  calculated. 


3.4  Wind  Speed 

The  wind  source  levels  are  computed  in  RANDI  3.1  with  data  extracted  from  the  HWS  data 
base  (Naval  Oceanographic  Office  1989).  The  HWS  data  base  contains  temporal  surface  statistics 
taken  from  marine  ships.  Included  are  the  number  of  reports  of  wind  speed  and  direction,  the  mean 
and  standard  deviation  of  the  wind  speed,  the  percentage  of  wind  direction  reports  for  calm  and 
variable  sea  states,  and  the  mean  wind  speed  and  its  wind  direction  for  45°  sectors.  The  HWS  data 
base  covers  the  Northern  Hemisphere  and  the  Southern  Hemisphere  to  65°  S.  The  resolution  is  a 
1°  grid  and  the  temporal  resolution  is  monthly.  This  data  base  is  unclassified. 


4.0  PROPAGATION  OF  SHIPPING  NOISE 

Once  a  source  has  been  located  and  the  source  level  determined  it  is  necessary  to  propagate 
energy  to  the  receiver  array.  The  objective  is  to  calculate  the  complex  pressures  at  the  hydrophones 
of  a  receiver  array  (Sec.  6.0).  One  way  to  find  the  answer  is  to  obtain  numerical  solutions  to  the 
wave  equation,  which  is  a  second-order  differential  equation.  However,  through  approximations,  it 
is  possible  to  obtain  a  parabolic  differential  equation  which  is  first  order,  and  it  is  much  easier  to 
solve  numerically.  Two  PE  methods  will  be  described  that  lead  to  a  calculation  of  the  pressures  at 
the  receiver  array.  One  of  the  PE  models  is  based  on  a  split-step  algorithm  and  is  referred  to  as  the 
SSPE  or  NSPE.  The  other  PE  model  is  based  on  a  finite  element  code,  and  it  is  referred  to  as  FEPE. 

The  environmental  data  (Sec.  3.0)  required  for  the  propagation  models  are  bathymetry,  sound 
speed  in  the  water,  sediment,  and  bottom;  and  attenuation  and  density  for  the  sediment  and  bottom 
layers.  The  bottom  characteristic  data  base  includes  geoacoustic  parameters  that  must  be  converted 
to  FEPE  or  SSPE  model  sediment  and  bottom  inputs. 


4.1  The  Parabolic  Approximation 

Underwater  acoustic  propagation  for  a  monochromatic  source  is  given  by  the  Helmholtz  equation, 
also  known  as  the  reduced  wave  equation.  In  cylindrical  coordinates  it  can  be  expressed  as 

(r,  (j),  z)p  =  0  ,  (7) 

where  p  is  the  acoustic  pressure,  ko  is  the  wave  number  given  by  (o/Co,  with  to  being  the  angular 
frequency  of  the  source,  and  Cq  is  the  sound  velocity  at  some  reference  depth.  The  index  of 
refraction  n  =  c/c^,  where  c  is  a  function  of  the  horizontal  range  r,  the  depth  z,  and  the  azimuthal 
angle  (j). 

Generally,  azimuthal  variations  of  the  sound  speed  are  small  and  can  be  neglected.  In  a  waveguide 
such  as  the  ocean,  cylindrical  spreading  occurs,  and  it  can  be  factored  out  with  the  substitution  of 


p(r,  z)  =  u(r,  z)e''^o^jyfr 


(8) 
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into  the  Helmholtz  equation.  The  phase  factor  &xp{ikor)  has  also  been  factored  out,  leaving  an 
envelope  function  u{r,  z)  that  obeys  the  parabolic  approximation  to  the  wave  equation 

likoduldr  +  d'^uldz^  +  z)  -  1)m  =  0  .  (9) 

As  can  be  seen  from  the  above  equation,  the  parabolic  approximation  is  a  first-order  differential 
equation  in  the  horizontal  range  r,  making  it  possible  to  numerically  solve  it  by  marching  an  initial 
field  out  in  range.  Two  approaches  to  solving  the  parabolic  equation  are  the  split-step  algorithm 
SSPE  and  a  finite  element  technique  leading  to  the  FEPE. 


4.2  Split-Step  PE 

The  split-step  algorithm  was  the  first  numerical  method  applied  to  the  solution  of  the  parabolic 
equation  by  Tappert  (1974),  and  it  involves  a  marching-type  Fourier  transform  solution.  Since  one 
cannot  numerically  deal  with  an  infinite  transform,  a  false  bottom  is  added  where  the  envelope  u  is 
attenuated.  The  SSPE  used  in  RANDI  3.1  is  a  wide-angle  (up  to  40°  half-beamwidth)  extension  of 
Tappert’s  original  solution  due  to  Thomson-Chapman  (1983),  and  it  has  been  adopted  as  the  NSPE. 

4.2.1  Characteristics  of  SSPE 

In  addition  to  its  relatively  wide-angle  capability,  the  SSPE  is  an  efficient  algorithm  to  compute. 
This  efficiency  comes  about  through  the  use  of  the  Fourier  transform  in  marching  the  solution  in 
range.  In  addition,  the  SSPE  uses  the  maximum  range  step  allowed  before  significant  errors  are 
introduced. 

The  SSPE  requires  small  range  steps  when  large  propagation  angles  are  involved.  In  addition, 
a  uniform  depth  grid  is  required  with  the  number  of  points  in  that  grid  confined  to  a  power  of  2 
for  speed  in  calculating  the  Fourier  transform.  Finally,  density  variations  or  discontinuities  are  a 
special  problem,  and  a  “smearing  out”  of  such  changes  is  necessary.  For  a  density  discontinuity 
from  pi  to  p2  at  the  depth  z  =  Zf,,  the  discontinuity  is  replaced  with  the  hyperbolic  tangent  smoother 

p(z)  =  pi  +  5(p2  -  pi){l  +  tanh[(z  -  zb)IL\}  ,  (10) 

where  L  defines  the  length  of  the  transition  region  over  which  the  discontinuity  is  smoothed, 
usually  chosen  such  that  kpL  =  2. 

4.2.2  SSPE  Inputs  from  Extracted  Bottom  Parameters 

An  advantage  of  the  NSPE  model  is  that  it  allows  the  user  to  input  the  BLUG  parameters 
(Sec.  3.0)  as  they  are  extracted.  The  values  are  converted  to  the  necessary  bottom  inputs,  such  as 
sediment  and  basement  sound  speed,  density,  and  attenuation.  The  only  input  that  must  be  derived 
is  the  sediment  thickness  Ds,  which  is  inferred  from  the  two-way  travel  time  T2.  The  relationship 
between  T2  and  Ds  is 


(1/2)72=1^  Cb{z)^dz, 


(11) 


where  in  BLUG,  the  sound  propagation  speed  in  the  bottom  Cb{z)  is  assumed  to  be  of  the  form 


cb  (z)  =  co{(i  +  m  +  2gz/(co(i  +  m''"-  p} , 


(12) 
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where  cq  is  the  initial  sediment  sound  speed,  g  is  the  sound  speed  gradient  at  the  water-sediment 
interface,  and  p  is  the  sediment  sound  speed  profile  curvature  factor. 

The  closed  form  solution  to  the  above  integral  is 

T2  =  2[(c£)/co)  +  P/m(co/co)  -  l]/[go(l  +  P)]  •  (13) 

This  equation  is  then  solved  numerically  for  C£>.  The  value  of  Dg  is  then  determined  from  the 
equation  for  Cb{z).  In  RANDI  3.1,  the  value  of  Ds  is  not  allowed  to  be  greater  than  2000  m. 


4.3  Finite  Element  PE 

Recently,  an  FEPE  model  (Collins  1988,  1989)  based  on  a  Fade  series  expansion  of  the 
depth-dependent  part  of  the  parabolic  equation  was  developed  in  an  effort  to  handle  propagation 
angles  greater  than  the  40°  half-beamwidth  capability  of  the  SSPE.  In  the  FEPE  model,  the  parabolic 
equation  is  solved  using  Galerkin’s  finite  elements  method  for  depth  discretization  and  the  Crank- 
Nicholson  method  for  integration  in  range.  This  code  is  used  in  RANDI  3.1,  and  it  contains  an 
efficient  tridiagonal  system  solver  designed  to  minimize  computation  time. 

4.3.1  Characteristics  of  FEPE 

The  tridiagonal  solver  is  more  efficient  than  the  Gaussian  elimination  method  for  range-dependent 
bathymetries.  Because  the  FEPE  model  does  not  require  a  uniform  depth  grid,  density  and  sound 
speed  discontinuities  can  be  better  approximated  without  a  substantial  increase  in  the  number  of 
depth  grid  points  as  in  the  SSPE.  However,  the  current  FEPE  code  does  not  handle  variable  depth 
grids.  The  Fade  series  approximation  to  the  depth-dependent  part  of  the  parabolic  equation  enables 
the  FEPE  to  propagate  at  angles  up  to  90°  half-beamwidth. 

Whereas  the  SSPE  can  be  regarded  as  automated  because  of  its  own  selection  of  range  step  and 
depth  transform  size,  the  FEPE  requires  that  the  user  input  both  range  and  depth  steps.  As  a  general 
rule,  the  maximum  range  increment  used  should  be  a  half-wavelength,  and  the  maximum  depth 
increment  used  should  be  a  quarter-wavelength.  The  increments  need  to  be  small  enough  to  adequately 
model  variations  in  the  environment.  The  FEPE  has  a  choice  of  starter  types  including  the  Gaussian 
starter,  the  Green’s  starter,  the  homogeneous  mode  starter,  and  the  image  starter. 

4.3.2  FEPE  Inputs  from  Extracted  Bottom  Parameters 

It  is  necessary  that  the  data  extracted  from  the  LFBL  data  base  (Sec.  3.0)  be  converted  to  values 
that  can  be  input  into  the  model.  As  in  the  SSPE  preparation,  Ds  is  calculated  using  the  two-way 
travel  time  and  the  expression  for  Cb{z)  given  in  Sec.  4.2.2.  The  sediment  sound  speed  profile  is 
calculated  using  the  equation  for  c^(z).  The  compressional  wave  velocity  in  the  basement  is  simply 
Cb{Ds),  where  the  sediment  density  is  entered  as  extracted  from  the  data  base.  The  basement  density 
is  calculated  using  the  relationship 

Pfc  =  p,(l+R)/(1-R),  (14) 

where  p^  is  the  sediment  density  (from  BLUG)  and  R  is  the  BLUG  basement  reflection  coefficient. 
The  attenuation  as  a  function  of  depth  is  calculated  from 


(2)  =  (//1000)(ct  iz)lf){ko  +  kQz)  , 


(15) 
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where  /  is  the  frequency  in  hertz,  (z)  is  the  sediment  sound  speed  at  depth  z,  ko  is  the  surface 
sediment  attenuation  factor  in  dB/m/kHz,  and  ko  is  the  attenuation  gradient  in  dB/m/kHz/m. 

If  a  Hamilton  type  bottom  is  used  in  the  model,  the  basement  is  assumed  to  be  basalt.  The 
compressional  wave  velocity  in  the  basement  is  estimated  to  be  5300  m/s.  The  density  in  the  sediment 
is  calculated  using  the  expression 

Pd  =  1.2  +  [(c,  -  cfe(zi))/1000]^/4  ,  (16) 

where 

Cs  =  {[cb(zi)][cb(Ds)]}^^2^  (17) 

The  density  in  the  basement  is  set  to  2.7  gm/cm^. 

4.3.3  FEPE  Example 

An  example  of  the  propagation  of  sound  determined  by  the  FEPE  propagation  model  is  shown 
in  Fig.  7.  The  sound  is  generated  by  five  ships  on  a  slope  in  the  Arabian  Sea.  The  ships  are  assumed 


Fig.  7  —  Example  of  downslope  noise  enhancement  as  calculated  by  the  FEPE  model  for  five  ships  on  a 

track  in  the  Arabian  Sea 
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to  have  equal  source  levels  and  are  spaced  15  km  apart.  The  location  of  each  source  is  seen  on  the 
left  side  of  the  figure.  The  legend  at  the  bottom  right  of  the  figure  indicates  the  signal  strength  in 
decibels.  In  this  example,  the  effects  of  downslope  enhancement  are  clearly  evident.  The  high-angle 
energy  from  the  ships  has  been  converted  to  long-range,  refracted-refracted  paths  cycling  at  a  depth 
of  1750  m. 


4.4  FEPE  or  SSPE? 

In  considering  the  choice  of  FEPE  or  SSPE,  the  user  must  evaluate  the  environment  being 
modeled.  The  FEPE  model  is  best  for  moderate  range  dependence,  it  handles  the  bottom  better,  and 
is  accurate  up  to  90°  half-beamwidth.  The  SSPE  model  is  best  for  weak  range  dependence,  is 
accurate  up  to  40°  half-beamwidth,  but  smoothes  density  discontinuities. 


5.0  WIND  NOISE 

A  major  source  of  ambient  noise  in  the  ocean  is  due  to  breaking  water  waves  (Felizardo  and 
Melville  1995).  Since  water  waves  are  primarily  generated  by  the  wind,  this  is  referred  to  as  wind 
noise.  Wind  noise  in  the  vicinity  of  the  receiver  array  is  computed  in  RANDI  3.1  using  a  model 
due  to  Kuperman  and  Ingenito  (1980)  which  was  extended  by  Hamson  (1985)  to  include  a  more 
general  treatment  of  source  directionality  functions.  A  brief  description  of  the  wind  model  follows. 


5.1  Kuperman-Ingenito  Wind  Noise  Model 

The  Kuperman-Ingenito  wind  model  is  range  independent  with  three  horizontal  layers.  The  top 
layer  is  the  water  column,  the  middle  layer  is  composed  of  sediments,  and  the  bottom  layer  is  the 
rest  of  the  ocean  bottom  extending  to  infinity.  In  both  the  water  and  sediment  layers  the  sound 
velocity  can  vary  with  depth.  In  the  bottom  two  layers  constant  values  of  density  and  attenuation 
are  assumed.  Since  the  wind  model  is  range  independent,  the  environmental  input  information  is 
extracted  only  at  the  receiver  location. 

An  infinite  layer  of  sources  is  assumed  to  exist  just  below  the  ocean  surface.  A  complete 
solution  to  the  wave  equation  is  found  which  includes  a  discrete  spectrum  (normal  modes)  and  a 
continuous  spectrum.  The  discrete  spectrum  results  from  waves  which  are  totally  reflected  at  the 
bottom  and  become  trapped  in  the  waveguide  created  by  the  layering.  The  spectrum  is  discrete 
since  energy  is  found  only  for  those  frequencies  where  constructive  interference  occurs  in  the 
waveguide.  The  continuous  spectrum  occurs  since  the  receiver  is  close  to  the  source  where  waves 
reflect  from  the  bottom  at  less  than  the  critical  angle  for  total  reflection.  Using  the  Kuperman- 
Ingenito  wind  model,  Hamson  (1985)  shows  that  continuous  modes  can  be  necessary  to  accurately 
model  measured  data  in  shallow  water,  especially  if  the  bottom  is  soft. 

Kuperman  and  Ingenito  (1980)  and  Hamson  (1985)  derive  expressions  to  compute  the  spatial 
correlation  function  of  the  noise  between  any  two  points.  The  discrete  mode  solution  is  obtained 
using  the  SACLANTCEN  Normal  Mode  Acoustic  Propagation  Model  (Jensen  and  Ferla  1979) 
and  the  continuous  field  solution  is  computed  using  the  Fast  Field  Program  (Kutschale  1973).  In 
RANDI  3.1,  the  wind  source  levels  are  determined  following  Wilson  (1983).  The  source  level  is  a 
function  of  the  frequency  and  wind  speed.  The  computed  wind  noise  is  obtained  as  a  complex 
spatial  correlation  matrix  between  all  the  hydrophones  in  the  receiver  array. 
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5.2  Wind  Noise  Bottom  Parameters 

The  value  of  the  sediment  thickness,  Ds,  is  computed  based  on  the  two-way  travel  time  in  the 
sediment  as  described  in  Sec.  4.2.2.  To  maintain  a  reasonable  number  of  modes  for  the  wind  noise 
calculation  in  RANDI  3.1,  Ds  is  not  allowed  to  be  greater  than  500  m. 

The  density  in  the  sediment  is  calculated  using  the  expression 

p,  =  1.2  +  [(c,  -  C6(zi))/1000]^/4  ,  (18) 

where  is  the  sediment  density,  is  the  sediment  sound  speed  at  depth  zi,  and 

Cs  =  {[cb(zi)][cb(Ds)]}^/^.  (19) 

The  remaining  bottom  parameters  refer  to  the  Hamilton  bottom  type  assuming  a  basalt  basement. 
Sediment  compressional  attenuation  is  0.38  dB/X,  sediment  density  is  2.7  gm/cc,  basement  com- 
pressional  attenuation  is  0.106  dB/X,  basement  compressional  speed  is  5300  m/s,  basement  shear 
attenuation  is  0.188  dB/X.  and  basement  shear  speed  is  2680  m/s. 


6.0  NOISE  PREDICTIONS  INDEPENDENT  OF  RECEIVER 

In  this  section,  the  prediction  of  ambient  noise  due  to  shipping  and  the  local  wind  will  be 
considered.  Three  methods  will  be  described  for  calculating  shipping  noise.  In  addition,  examples 
will  be  presented  of  noise  predictions  which  are  independent  of  the  receiver  array. 

In  the  RANDI  3.1  model,  the  shipping  complex  pressures  are  computed  by  executing  one  of 
the  PE  propagation  models,  either  the  FEPE  model  or  the  SSPE  model  (Sec.  4.0).  To  begin,  the 
model  must  have  a  shipping  field  defined  by  ship  locations  and  strengths.  The  environmental 
shipping  density  data  bases  (HITS  3.0  and  3.1)  are  used  to  retrieve  the  shipping  densities,  and  these 
densities  are  converted  to  discrete  ship  locations.  Source  levels  are  assigned  depending  on  the  class 
of  ship  (Sec.  2.0).  The  environmental  data  required  for  the  propagation  models  are  bathymetry; 
sound  speed  in  the  water,  sediment,  and  bottom;  and  attenuation  and  density  for  the  sediment  and 
bottom  layers  (Sec.  3.0).  The  bottom  characteristic  data  base  includes  geoacoustic  parameters  that 
must  be  converted  to  FEPE  or  SSPE  model  sediment  and  bottom  inputs.  Bottom  parameters  derived 
from  Hamilton  (1980)  (ex.  in  Sec.  5.2)  can  also  be  included  in  the  sediment  and  bottom  inputs 
(Secs.  3.0  and  4.0). 

In  considering  the  choice  of  FEPE  or  SSPE  models,  the  user  must  evaluate  the  environment 
being  modeled  and  the  location  of  the  sources  with  respect  to  the  receiver.  For  example,  are  the 
sources  near  or  distant?  The  SSPE  model  can  produce  accurate  results  and  be  computationally 
efficient  where  the  propagation  is  narrow  angle  (with  respect  to  the  horizontal)  and  weakly  range 
dependent.  However,  discontinuities  at  the  water-bottom  interface  create  problems  for  the  SSPE 
model.  For  complicated  bottom  interactions,  a  small  computational  grid  is  required  to  obtain  accu¬ 
rate  results,  and  this  can  severely  reduce  the  computational  efficiency  of  the  algorithm.  While  the 
SSPE  model  is  accurate  for  up  to  40°  half-beamwidth,  the  FEPE  model  is  accurate  up  to  90°  half- 
beamwidth.  Therefore,  the  FEPE  model  is  much  more  accurate  than  the  SSPE  model  in  predicting 
the  noise  of  nearby  sources.  The  FEPE  is  also  the  more  accurate  of  the  two  models  for  range- 
dependent  environments  such  as  found  in  shallow  water.  Also,  in  this  situation,  the  SSPE  offers  no 
computational  advantage  (Jensen  et  al.  1994;  Etter  1991). 
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6.1  Methods  for  Calculating  Shipping  Noise 

In  RANDI  3.1,  three  different  methods  are  provided  for  computing  shipping  noise.  These 
methods  are  referred  to  as  the  rigorous  method,  the  radial  accumulation  method,  or  the  high- 
resolution  radial  accumulation  method.  Each  method  is  described  below. 

6.1.1  Rigorous  Method 

The  simplest  and  most  accurate  method  of  computing  shipping  noise,  which  is  referred  to  as 
the  rigorous  method,  entails  running  one  of  the  PE  models  from  each  ship  to  the  receiver.  The 
complex  pressures  are  saved  over  a  portion  of  the  range-depth  mesh  which  encompasses  the  receiver. 
Obviously,  this  method  involves  as  many  PE  runs  as  there  are  ships,  and  can  be  very  time-consuming 
for  long  ranges  and  a  large  number  of  ships.  The  number  of  bathymetric  and  sound  speed  profile 
extractions  needed  is  also  proportional  to  the  number  of  ships. 

6.1.2  Radial  Accumulation  Method 

One  way  to  reduce  the  number  of  environmental  extractions  and  propagation  calculations  required 
is  to  discretize  the  azimuth  into  a  number  of  “pie-slice”  sectors  which  are  referred  to  as  radial 
sectors.  This  is  called  the  radial  accumulation  method,  and  it  works  only  with  the  FEPE  model.  As 
shown  in  Fig.  8,  ships  are  now  assigned  to  particular  radial  sectors  depending  on  their  locations. 

The  ship  positions  in  a  sector  are  projected  onto  a  radial  selected  for  the  sector,  and  the  FEPE 
is  run  from  the  farthest  ship  to  the  receiver,  accumulating  other  ships  as  it  propagates  along  the 
selected  radial  toward  the  receiver.  For  the  particular  radial,  the  total  pressure  is  given  by 

Pj^Uje^^^lypTc,  (20) 

where  Vc  is  the  horizontal  range  from  the  farthest  ship  in  that  radial  sector  to  the  receiver  center  and 

ships 

Ut=  i  LsnUn>JrJ^e-^Krc-r„)  (21) 

n  =  1 


TL  CALCULATION  FROM 
SHIP  TO  RECEIVER 


Fig.  8  —  Diagram  of  radial  accumulation  method  for  propagation  of  ships 
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is  the  total  envelope  function  obtained  by  accumulating  the  ships  with  source  levels  L^n  as  the 
FEPE  propagates.  This  approximation  to  the  total  pressure  for  a  radial  is  valid  only  in  a  local  region 
surrounding  the  receiver. 

The  above  algorithm  is  repeated  for  each  radial  sector,  and  it  can  be  seen  that  the  number  of 
FEPE  runs  has  been  reduced  from  the  number  of  ships  to  the  number  of  radial  sectors  into  which 
the  azimuth  has  been  discretized,  saving  considerable  computing  time.  However,  there  is  a  sacrifice 
in  the  accuracy  of  the  results,  which  decreases  with  increasing  sector  size.  When  dealing  with 
thousands  of  ships,  the  rigorous  method  can  become  impractical  because  of  the  excessive  compu¬ 
tational  time.  By  comparison,  reasonable  results  can  often  be  obtained  with  much  less  computational 
time  by  using  the  radial  accumulation  method  with  small  sectors. 

6.1.3  High-Resolution  Radial  Accumulation  Method 

To  regain  some  of  the  individual  ship  information  lost  by  the  radial  accumulation  method, 
RANDl  3.1  also  includes  an  option  called  the  high-resolution  radial  accumulation  method.  The  ships 
are  accumulated  along  a  radial  in  a  sector  just  as  described  in  the  previous  section  using  the  FEPE 
model.  Weighting  coefficients  are  then  estimated  for  each  ship  in  the  sector,  and  these  coefficients 
are  used  to  estimate  the  complex  pressures  of  the  corresponding  ships.  For  a  particular  sector,  the 
weighting  coefficient  W5  for  a  ship  s  is  given  by 

20  \ogWs  =  SLs-TLs,  (22) 

where  SLs  is  the  signal  level  of  ship  s  at  the  receiver  and  TLs  is  its  transmission  loss  to  the  receiver 
center  as  approximated  by  10  log  r^,  where  is  the  range  from  the  ship  to  the  center  of  the 
receiver.  The  weighting  coefficients  are  normalized  by  dividing  by  the  sum  of  all  the  weights  in 
the  sector,  i.e.,  by  the  factor 


ships 
5=  1 


(23) 


With  these  weighting  coefficients,  each  ship  in  a  radial  sector  is  assigned  the  complex  pressure 

Ps  —  Pradial  >  (24) 

where  Pradial  is  the  total  complex  pressure  for  the  radial  sector  determined  by  the  radial  accumulation 
method.  In  Sec.  7.1.3,  where  the  interest  is  beamforming,  an  adjustment  in  phase  based  on  a  ship’s 
bearing  is  applied  to  the  pressure  at  the  receiver  (hydrophones)  for  each  ship. 

Under  most  circumstances,  the  high-resolution  radial  accumulation  method  would  be  chosen 
over  the  radial  accumulation  method,  since  the  computational  time  required  for  the  additional 
adjustments  in  complex  pressure  and  phase  is  negligible.  However,  the  high-resolution  adjustments 
are  only  approximations,  and  although  tests  have  shown  that  the  approximations  are  generally  good, 
for  some  situations  the  results  may  not  be  sufficiently  accurate,  or  may  even  be  misleading. 


6.2  Shipping  and  Wind  Noise 

Depending  upon  which  shipping  noise  calculation  option  is  chosen,  the  propagation  model  is 
then  run  from  the  ship  (or  maximum  range  of  the  radial)  to  the  receiver,  ignoring  ships  or  radials 
blocked  by  land,  and  saving  and  accumulating  complex  pressures  at  the  receiver.  The  contribution 
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due  to  the  local  wind,  if  desired,  is  determined  for  a  range-independent  environment  (Sec.  5.0)  and 
added  to  the  shipping  noise. 

The  horizontal  shipping  noise  field  in  the  East  China  Sea  near  the  Korea  Strait,  as  calculated 
by  the  three  propagation  methods,  are  shown  in  Figs.  9,  10,  and  11  for  a  frequency  of  30  Hz.  The 
maximum  ship  range  is  3000  km,  but  the  number  of  ships  has  been  randomly  thinned  to  decrease 
computational  time.  The  receiver  is  located  at  31°  N,  127°  E.  These  “spiked”  plots  represent  the 
horizontal  noise  field  without  the  array  response,  i.e.,  independent  of  the  receiver  particulars  such 
as  number  of  hydrophones  or  hydrophone  spacing. 

Figure  9  gives  the  noise  field  as  calculated  by  the  rigorous  method.  Each  spike  represents  the 
ship  noise  contribution  with  1°  bearing  resolution.  Nulls  in  the  field  are  areas  of  negligible  or 
nonexistent  shipping,  blockage  by  land,  or  a  lack  of  environmental  data.  It  is  possible  to  provide 
the  model  with  environmental  data  when  they  are  missing  in  the  included  data  bases.  However, 
implementation  of  the  data  must  be  performed  on  a  case-by-case  basis  for  a  given  site  and  data 
base. 

Figure  10  shows  the  horizontal  shipping  noise  field  as  determined  by  the  radial  accumulation 
method.  The  resolution  is  equal  to  the  radial  sector  size  which  is  5°,  corresponding  to  72  equally 
spaced  radials.  By  comparison  with  Fig.  9  it  is  seen  that  the  general  character  of  the  noise  field  is 
preserved.  However,  a  compromise  in  resolution  is  evident.  At  a  bearing  of  about  45°,  two  spikes 
are  seen  close  together  in  Fig.  9.  These  two  spikes  become  one  stronger  spike  in  Fig.  10,  since  it 
is  not  possible  to  distinguish  between  bearings  in  a  common  sector. 

Figure  1 1  shows  the  comparable  results  obtained  using  the  high-resolution  radial  accumulation 
method.  Noise  levels  have  been  assigned  to  the  individual  ships  using  weighting  coefficients  as 
described  in  Sec.  6.1.3.  The  ships  have  been  assigned  their  original  bearings. 

It  is  interesting  to  note  that  184  ships  contributed  to  the  noise  field  when  using  the  rigorous 
method,  whereas  in  the  radial  methods,  there  were  212  ships  contributing  to  the  noise.  This  happened 
because  in  using  the  rigorous  method  the  propagation  path  is  checked  between  each  ship  and  the 
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Fig.  9  —  Horizontal  shipping  noise  field  for  30  Hz  in 
the  East  China  Sea  calculated  by  the  rigorous  propagation 
method 


Fig.  10  —  Horizontal  shipping  noise  field  for  30  Hz  in 
the  East  China  Sea  calculated  by  the  radial  accumulation 
propagation  method 


RAND  I  3.1  Physics  Description 


19 


receiver.  A  number  of  ships  were  obviously  found  to  be  blocked  from  the  receiver.  In  using  the 
radial  methods,  ships  in  a  sector  are  only  ignored  if  blockage  is  found  along  the  propagation  radial 
specified  for  the  given  sector. 

Figure  12  depicts  the  horizontal  noise  field  of  Fig.  9  with  noise  from  a  10-kt  wind  included. 
The  omnidirectional  level  of  50.36  dB,  as  calculated  by  the  Kuperman-Ingenito  model,  is  divided 
among  360  1°  sectors  resulting  in  a  minimum  noise  level  of  24.80  dB  for  all  azimuthal  angles. 

In  RANDI  3.1,  the  bathymetry  and  sound  speed  profiles  extracted  for  use  in  the  propagation 
model  are  saved.  In  addition,  the  pressure  field  for  a  user-selected  ship  or  radial  is  saved.  For 
example.  Fig.  13  contains  the  environmental  information  from  the  5°  radial  sector  centered  on 
222.5°  used  for  the  calculation  of  the  noise  by  the  three  methods.  Five  ships  contribute  to  the  field 
in  this  sector.  The  corresponding  pressure  field  calculated  by  the  FEPE  is  shown  in  Fig.  14. 
Figure  15  depicts  the  FEPE  signal  level  (5L  -  TL)  versus  range  from  the  farthest  ship  to  the  receiver. 
The  first  four  ships  are  located  between  0  and  37  km.  The  fifth  ship  is  located  at  approximately 
458  km,  and  it  produces  a  noticeable  jump  in  the  signal  level  where  its  source  level  is  added  to  the 
field,  as  can  be  observed  in  Figs.  14  and  15. 


6.3  Dead-Reckoning 

The  RANDI  3.1  model  has  an  option  to  dead-reckon  the  ships  and  compute  the  shipping  noise 
contributions  for  each  dead-reckoning  period.  The  ships  are  directed  to  move  along  shipping  lanes. 
This  is  a  direction  assumed  to  be  perpendicular  to  the  gradient  of  the  shipping  density,  which  is 
determined  by  a  quadratic  fit  to  twelve  1°  grid  cells  surrounding  the  ship.  Shipping  lanes  are 
assumed  to  be  where  the  greatest  concentration  of  ships  are  found.  Once  all  the  time  periods 
are  modeled,  an  average  horizontal  noise  field  is  computed  by  averaging  each  spatially  smoothed 
noise  field  calculated  for  each  dead-reckoning  period.  An  example  is  shown  in  Fig.  16.  The  ships 
surrounding  the  Korea  Strait  are  dead-reckoned  twice  over  1-h  time  periods.  The  high-resolution 
radial  accumulation  method  is  used  to  calculate  the  shipping  noise  field.  The  local  wind  noise 
contribution  (50.36  dB)  is  added  before  smoothing  and  time-averaging  is  performed. 
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Fig.  11  — Horizontal  shipping  noise  field  for  30  Hz  in  Fig.  12  —  Horizontal  noise  field  due  to  ships  and  a 
the  East  China  Sea  calculated  by  the  high-resolution  10-kt  wind 

radial  accumulation  propagation  method 
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Fig.  13  —  Bathymetry  and  sound  speed  profiles  for  the  radial  sector  centered  on  222.5° 


7.0  NOISE  PREDICTIONS  DEPENDENT  OF  RECEIVER 

The  RANDI  3.1  model  produces  noise  predictions  for  one  of  several  receiver  array  types: 
horizontal,  vertical,  tilted  linear,  or  volumetric.  Linear  arrays  are  assumed  to  be  undeformed  and 
absolute  element  locations  are  specified  via  an  array  center  depth,  element  spacing,  array  heading, 
and  tilt  angle.  Hydrophone  locations  in  volumetric  arrays,  of  which  the  linear  array  is  a  subset,  are 
input  in  a  rigorous  manner  using  a  local  Cartesian  coordinate  system  at  the  receiver  referenced  to 
the  array  center  depth. 

The  focus  in  this  section  will  be  on  predictions  of  noise  that  depend  upon  the  array  of  sensors 
(hydrophones)  that  are  used  to  measure  the  noise.  It  should  be  noted  that  measurement  arrays  do 
not  directly  measure  a  noise  property  such  as  directionality.  Rather,  they  sample  the  acoustic 
pressure  at  individual  hydrophone  locations.  Beamforming  is  then  used  to  estimate  the  direction¬ 
ality  of  the  ambient  noise  field  at  a  receiver  site  as  a  function  of  arrival  angle  about  the  receiver 
location. 

One  of  the  most  powerful  aspects  of  the  RANDI  3.1  noise  model  is  that  it  can  be  used  to 
accurately  model  the  coherent  and  incoherent  beam  response  of  a  measurement  array  to  the  shipping 
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and  wind  noise  fields  at  a  receiver  location.  Additionally,  it  is  possible  to  predict  the  measurement 
array  responses  to  system  and  flow  noises. 


7.1  Shipping  Noise:  Complex  Pressures 

It  is  necessary  to  estimate  the  complex  pressures  about  an  array  before  it  is  possible  to  beamform 
the  results.  The  shipping  noise  is  computed  using  either  the  FEPE  or  the  SSPE  model  (Sec.  4.0) 
with  one  of  three  methods:  the  rigorous  method,  the  radial  accumulation  method,  or  the  high- 
resolution  radial  accumulation  method  (Sec.  6.0).  Regardless  of  which  method  is  employed,  the 
values  of  the  envelope  function  u(r,  z),  which  obeys  the  parabolic  approximation  to  the  wave 
equation  (Sec.  4.1),  are  saved  for  range-depth  gridpoints  in  the  vicinity  of  the  receiver  array  of 
hydrophones.  In  general,  the  position  of  these  gridpoints  will  not  coincide  with  the  positions  of  the 
array  hydrophones,  and  some  type  of  interpolation  will  be  necessary  to  obtain  the  complex  pressures 
at  each  hydrophone. 

7.1.1  Parameters  for  a  Tilted  Array 

The  calculation  of  the  complex  pressure  at  a  hydrophone  n  makes  use  of  the  fact  that  the 
envelope  function  u{r,  z)  calculated  by  the  PE  is  a  slowly  varying  function  of  range.  As  a  consequence 
of  this  assumption,  a  linear  interpolation  between  the  gridpoints  is  a  good  approximation.  Toward 
this  aim,  both  the  range  and  depth  coordinates  of  each  hydrophone  are  calculated. 

To  be  general,  it  will  be  assumed  that  the  array  is  tilted  at  some  angle.  The  array  tilt  is  specified 
with  the  use  of  a  local  coordinate  system  with  hydrophone  1  at  its  origin.  The  array  elevation 
angle  0  is  measured  from  the  horizontal  with  the  convention  that  it  is  positive  in  the  direction  of 
increasing  z,  i.e.,  depth,  and  its  heading  (j)  is  referenced  from  north,  being  positive  in  the  clockwise 
direction.  Given  the  separation  d  between  hydrophones,  the  depth  of  hydrophone  n  is  given  by 

2^  =  ZQ  +  (rt  -  l)d  sin  0  ,  (25) 

so  that  hydrophone  1  is  always  located  at  the  origin  whose  depth  is  zq,  and  n  runs  from  1  to 
N  hydrophones. 

The  range  used  to  calculate  u(r,  z)  is  the  distance  from  each  individual  ship  to  the  receiver  for 
the  rigorous  method,  and  the  distance  from  the  farthest  ship  on  the  radial  to  the  receiver  for  the 
radial  accumulation  methods.  Because  the  ranges  involved  in  the  PE  calculations  are  rather  large, 
a  plane  wave  approximation  is  used  in  the  estimation  of  the  range  to  each  hydrophone  n 

=  ro  -  (n  -  l)d  cos  0  cos  ((^s  -  ^)  ?  (26) 

where  ro  is  the  range  separation  between  the  origin  of  the  local  array  coordinate  system  and  each 
ship  and  (j)^  is  the  ship’s  bearing  relative  to  north  as  seen  in  the  array  coordinate  system.  For  the 
radial  accumulation  method,  ro  is  the  range  to  the  farthest  ship  in  each  sector  and  (^s  is  the  bearing 
of  the  sector  radial  with  respect  to  north,  again  as  observed  from  the  array. 

7.1.2  Interpolation  of  the  Complex  Pressures 

The  range-depth  gridpoints  at  which  the  PE  runs  produce  values  for  the  envelope  function 
i/(r,  z)  and  are  separated  in  range  by  a  maximum  of  a  half- wavelength  (SSPE  and  FEPE).  The 
separation  in  depth  is  a  maximum  of  a  quarter-wavelength  for  the  FEPE  model,  or  a  fraction  of  a 
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wavelength  (determined  by  the  size  of  the  Fourier  transform)  for  the  SSPE  model.  The  complex 
pressures  at  each  hydrophone  are  obtained  by  linearly  interpolating  the  values  of  the  envelope 
u{r,  z)  over  each  hydrophone’s  range  and  depth  and  then  multiplying  the  interpolated  value  of  u  by 
a  factor  which  incorporates  the  phase  due  the  horizontal  range  and  cylindrical  spreading.  For  a 
hydrophone  of  index  n,  the  interpolated  complex  pressure  is  given  by 

^n)  =  «('•«>  /  \/^  ’  (27) 

where  is  given  by  the  equation  in  the  preceding  section,  and  Zn  is  the  azimuthally  independent 
depth  of  hydrophone  n. 

7,1.3  Phase  Estimation  for  the  High-Resolution  Radial  Accumulation  Method 

Weighting  coefficients  were  obtained  in  Sec.  6.1.3  to  assign  pressures  to  individual  ships  when 
the  high-resolution  radial  accumulation  method  is  used.  Before  it  is  possible  to  do  beamforming, 
it  is  also  necessary  to  apply  an  adjustment  in  phase  to  each  ship  pressure  in  a  radial  sector.  Let  all 
bearings  be  referred  to  north  and  be  positive  in  the  clockwise  direction.  For  a  radial  whose  median 
has  the  bearing  the  phase  delay  between  hydrophones  for  a  plane  wave  reaching  a  horizontal 
array  with  heading  (t)^^  is  equal  to 


^am  —  ho  d  cos 

where  is  the  wave  number  and  d  is  the  separation  between  hydrophones.  The  actual  bearing  of 
ship  5,  however,  is  (j)^,  so  that  a  plane  wave  emanating  from  that  bearing  would  suffer  a  phase  delay 
across  this  horizontal  array  equal  to 

aas  ^kod  cos  {(\>a  -  <1)5)  •  (29) 

For  a  horizontal  array,  then,  the  complex  pressure  assigned  to  ship  5  is  given  by 

Ps  =  WaPiadia\  ,  (30) 

where  Aa  =  ko  ^/[cos  ((t)^  -  (^s)  -  cos  ((j)^  - 

The  above  result  for  a  horizontal  array  can  easily  be  generalized  to  a  tilted  array,  giving  the 
following  expression  for  the  complex  pressure  assigned  to  ship  s 

Ps  =  Wa/7radial  ,  (31) 

where  6^  is  the  elevation  angle  of  the  array  with  reference  to  the  horizontal  plane. 


7,2  Shipping  Noise:  Noise  Directionality 

Propagation  algorithms  based  on  ray  theory  lead  to  straightforward  calculations  of  received 
energy  versus  elevation  angle,  the  PE-based  algorithms  do  not.  However,  this  information  can  be 
derived  from  the  PE  field  at  the  receiver  by  sampling  the  water  column  in  the  vertical  direction  and 
then  beamforming  to  obtain  the  vertical  arrival  structure  for  each  ship  or  radial  bearing,  depending 
on  whether  the  rigorous  or  radial  accumulation  method  is  chosen.  In  what  follows,  we  restrict 
ourselves  to  the  latter  method. 
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7.2 J  Cosine-Spaced  Beams 

Beams  are  steered  in  the  vertical  so  that  they  are  evenly  partitioned  in  cosine  space.  If  the  beam 
index  b  runs  from  1  to  the  number  of  beams  B,  each  beam  is  steered  so  that 

cos  =  1  -  2(6  -  1)/^  ,  (32) 

where  0^  is  the  steering  angle  as  measured  from  the  array  heading.  The  desirable  feature  of  cosine¬ 
spaced  beams  is  the  fact  that  they  subtend  equal  areas  on  the  unit  sphere.  This  area  is 

=  l  J  =  ^(COS0)  (33) 


or 


^azimuth  ~  4jr/B  ,  (^4) 

which  is  independent  of  the  beam  index  b. 

For  the  radial  accumulation  method,  the  azimuth  is  broken  up  into  Nr  radials  so  that  the  area 
Ar  subtended  by  the  cosine-spaced  beams  over  a  radial  is 

Ar  =  Aazimuth!^ r  “  ^“KlBNy  •  (35) 

7.2.2  Three-Dimensional  Directionality 

Assume  that  each  cosine-spaced  beam  is  located  at  the  center  of  each  Ar-  Then  the  beam 
intensity  (power  per  steradian)  over  an  area  characterized  by  the  beam  index  b  and  the  radial  index 
r  is 


Ibr=\pbr\^  (36) 

where  pbr  is  the  sum  of  the  delayed  complex  pressures  over  the  array  for  radial  r.  The  time  delay 
imposed  on  the  complex  pressures  depends  on  the  beam  index  b  such  that 

Axn,  b  =  nd  cos  0b/c ,  (37) 

where  n  is  the  hydrophone  index,  d  is  the  spacing  between  elements,  and  c  is  the  reference  sound 
speed. 

From  the  above  BNr  values  of  Ibr,  a  panoramic  mural  of  the  directionality  of  the  total  ambient 
noise  can  be  created.  An  example  of  the  noise  three-dimensional  directionality  for  the  example 
described  in  Sec.  6.2  is  given  in  the  color  plot  of  Fig.  17.  The  color  scale  appears  in  the  bottom 
right  corner  of  the  figure.  Each  vertical  striation  of  color  represents  the  noise  level  for  one  radial. 

Noise  directionalities  for  the  vertical  and  horizontal  can  easily  be  obtained  by  appropriate 
summations  of  the  BNr  values  of  Ibr-  For  example,  the  vertical  directionality  of  the  ambient  noise 
is  obtained  by  computing  the  noise  intensity  impinging  on  the  center  of  the  array  at  a  steering  angle 
0^  from  all  radials. 

Recalling  from  the  discussion  above  that  all  cosine-spaced  beams  subtend  an  equal  area,  the 
vertical  beam  intensity  (power  per  steradian)  is 
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4=2  Pbr  p  /^azimuth  =  2  hrl^r  ’  (38) 

r  I  r 

where  the  last  equality  can  be  verified  by  substitution  of  the  expression  for  I^r  from  the  previous 
section.  The  vertical  beam  intensity  corresponding  to  a  steering  angle  0^,  is  therefore  obtained  by 
averaging  over  the  number  of  radials.  An  example  of  the  vertical  directionality  is  shown  on  the 
right  side  of  Fig.  17.  For  this  example,  there  is  little  variability  in  the  vertical  direction  for  a 
frequency  of  30  Hz  since  the  water  depth  is  less  than  two  wavelengths  deep.  If  the  receiver  array 
is  located  farther  from  the  coast  at  23°  N,  150°  E,  the  three-dimensional  and  vertical  directionalities 
change  as  shown  in  Fig.  18.  The  increased  noise  level  near  horizontal  arrival  angles  in  the  vertical 
arrival  structure  reflects  the  increased  distance  between  the  coastal  areas  of  heavy  shipping  activity 
and  the  receiver  array.  The  omnidirectional  noise  level  is  correspondingly  decreased  from  89.7  dB 
to  71.7  dB. 

The  horizontal  directionality  of  the  total  ambient  noise  is  another  easily  obtainable  quantity 
from  the  Ihr  values.  To  calculate  the  noise  intensity  impinging  on  the  center  of  the  array  at  an 
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Fig.  17  —  Three-dimensional,  horizontal,  and  vertical  noise  directionalities  for  a  receiver  located  in  the  East 

China  Sea  near  the  Korea  Strait;  ocean  depth  99  m 
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Fig.  18  —  Three-dimensional,  horizontal,  and  vertical  noise  directionalities  for  a  receiver  located  in  deep 
water  outside  the  East  China  Sea;  ocean  depth  5647  m 


azimuthal  angle  ({);•,  it  is  necessary  to  sum  the  over  the  beam  index  b  for  each  radial.  The  area 
thus  subtended  is  AsUce  =  so  that 


Pbr^  ! ^ slice  =  2  ’ 

b 


(39) 


where  the  last  equality  can  again  be  verified  by  substitution  of  the  expression  for  Ibn  The  horizontal 
beam  intensity  (power  per  steradian)  over  a  radial  is  therefore  obtained  by  averaging  over  the 
number  of  beams. 


The  horizontal  directionality  is  displayed  in  terms  of  power  per  degree  slice  instead  of  beam 
intensity  (power  per  steradian).  No  matter  what  the  angular  width  of  the  radials,  one  can  always 
obtain  the  horizontal  beam  power  over  1°  slices  by  dividing  the  sum  of  powers  over  the  beams  by 
the  angular  width  of  a  radial.  Then 


Power  =  Pbr 

b 


2 


#^/360  =4jt/,/360. 


(40) 
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Examples  of  the  horizontal  directionality  computed  by  this  method  are  depicted  at  the  bottoms 
of  Figs.  17  and  18.  As  expected,  the  azimuthal  angles  between  approximately  180°  and  270° 
contain  the  highest  levels  of  shipping  noise,  coming  mainly  from  the  East  China  Sea. 


7.3  Spatial  Filtering 

To  perform  spatial  filtering,  RANDI  3.1  includes  Discrete  Fourier  Transform  (DFT)  and  Fast 
Fourier  Transform  (FFT)  beamformers,  each  of  which  is  applied  in  conjunction  with  one  of  several 
standard  array  shading  functions.  The  noise  pressures  can  be  summed  either  coherently  or  incoherently 
during  beamforming. 

7.3.1  Shading  Coefficients 

In  beamforming  noise,  it  is  often  desirable  to  employ  shading  coefficients.  Four  beamformer 
spatial  shading  options  are  available  in  RANDI  3.1.  They  are  unshaded  (uniform),  Hann,  Hamming, 
and  Blackwell.  The  shading  coefficients  for  these  four  options  are  computed  as  follows: 

For  an  unshaded  array. 


a„  =  1 . 


For  a  Hann  shaded  array. 


an 


=  0.5  -  0.5  cos 


2jt/i 

N  +  \‘ 


For  a  Hamming  shaded  array, 

in  =  0.54  -  0.46  cos  . 

N-\ 

For  a  Blackwell  shaded  array, 

^  2Kn  „  „„  4nn 

a„  =  0.42  -  0.50  cos - +  0.08  cos - . 

N+1  N+1 

The  coefficients  an  are  then  defined  to  be 


Un  =  UnlNorm , 


where 


(41) 

(42) 

(43) 

(44) 

(45) 


N 

Norm  ^  Un  (46) 

W  =  1 


is  a  normalization  constant. 

7.3.2  Discrete  Beamformers 

In  RANDI  3.1,  the  beam  response  due  to  shipping  noise  for  a  receiver  array  can  be  computed 
using  either  a  coherent  or  an  incoherent  beamformer.  The  DFT  coherent  beamformer  is  defined  by 
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2 


B 


Ships,  b 


N 

2  a 

n  =  l 


(47) 


where  Bships,b  is  the  beam  noise  power  due  to  shipping  on  beam  b,  Pships,n  is  the  total  complex 
pressure  at  hydrophone  n  due  to  all  ships,  a„  is  a  shading  coefficient,  and  is  the  time  delay 
for  hydrophone  n  used  to  form  beam  b. 

The  DFT  incoherent  beamforming  is  computed  by 


B 


Ships,  b 


Ships 

=  2 

s  =  1 


N 

'S.  a„P,,,txp(i27tfAx„l, 


2 


(48) 


where  Ps,n  is  the  complex  pressure  at  hydrophone  n  due  to  ship  5. 

For  a  line  array  with  phones  spaced  a  distance  d  apart,  the  time  delay  Axn,b  is  given  by 

Ax„^h  =  nd  cos  9i,/c  (49) 

and  the  beams  are  cosine  spaced  as  described  in  Sec.  7.2.1,  where  B  is  the  number  of  beams. 

An  FFT  beamformer  is  also  available  in  RANDI  3.1.  The  output  power  for  beam  m,  m  =  l, 
2, ... ,  N,  of  the  FFT  beamformer  is  given  by 


B 


Ships,  m  ' 


N 

2  ^ Ships,  n 

«  =  1 


n-\ 


exp 


-/2jt(«  -  l)(m  -  1)/Ar] 


(50) 


The  FFT  and  the  DFT  beamformers  will  produce  identical  results  for  the  special  case  when 
B  =  N  and  the  receiver  array  elements  are  half-wavelength  spaced.  At  frequencies  below  the  design 
frequency  of  the  receiver  array,  the  FFT  beamformer  forms  virtual  beams  corresponding  to  beams 
with  time  delays  larger  than  necessary  to  form  an  endfire  beam.  Berrou  and  Wagstaff  (1982)  have 
shown  that  virtual  beams  can  be  useful  in  tow  system  analysis.  In  general,  the  beams  for  the  FFT 
beamformer  point  in  directions  which  satisfy 

cos  =  ifdlf)  [1  -  2(m  -  1)/A^.  (51) 

Note  that  the  beam  directions  for  the  FFT  beamformer  depend  upon  frequency.  Also,  the  FFT 
algorithm  used  in  the  beamformer  requires  that  AT  be  a  power  of  2. 


7.4  Array  Response  to  Wind  Noise 

In  contrast  to  shipping  noise,  wind  noise  is  a  random  process.  Local  wind  noise  is  computed 
using  the  Kuperman-Ingenito  model  (Sec.  5.0),  which  is  range  independent.  An  infinite  layer  of 
sources  is  located  just  below  the  surface  and  the  bottom  is  layered.  Both  discrete  (normal  modes) 
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and  continuous  spectra  are  obtained.  The  wind  model  output  is  a  cross-spectral  density  matrix, 
which  is  used  in  beamforming. 

The  beam  outputs  due  to  wind  noise  are  determined  by 


N 


B 


Wind,  b 


=  1 


%  ^Wind^nm 


+  2  2 

n>m 


a„a^Re 


(52) 


where  Cwind,mn  is  the  cross-spectral  density  matrix  due  to  wind  noise  at  each  element  of  the 
receiver  array  and 


At/,  =  d  cos  0/,/c  (53) 

is  the  inner  hydrophone  time  delay  required  to  steer  a  beam  in  direction  6/,. 


7.5  Array  Response  to  Self  Noise 

In  RANDI  3.1,  there  is  an  option  to  compute  self  noise,  which  consists  of  system  noise  and 
flow  noise.  Both  of  these  components  to  the  total  noise  field  are  described  in  terms  of  cross-spectral 
density  matrices. 


7.5.1  System  Noise 

The  effect  of  system  noise  on  beamforming  is  included  in  RANDI  3.1  by  specifying  the  receiver 
array  cross-spectral  density  matrix  due  to  system  noise.  The  elements  in  this  cross-spectral  density 
matrix  are  assumed  to  be  of  the  form 

^ System, mn  —  ^System  sinc[jc»id/D] ,  (54) 

where  Lsystem  is  the  omnidirectional  system  noise  at  the  hydrophone  level,  d  is  hydrophone  spacing, 
and  D  is  a  user-specified  correlation  length.  This  approach  allows  the  user  great  flexibility  in 
characterizing  system  noise.  For  djD  »  1,  system  noise  is  effectively  uncorrelated  from  hydro¬ 
phone  to  hydrophone.  For  dfD  «  1,  system  noise  is  correlated  from  hydrophone  to  hydrophone 
with  the  amount  of  correlation  defined  by  the  sine  function.  Given  the  cross-spectral  density  matrix 
Csystem,mn  ^^d  the  desired  steering  direction  9/,,  the  beam  output  for  system  noise  is  given  by 

N 

2 

^ System,  b  ~  ^  System  ^  ^  ^m^5y5K/n,/nn^®^P^/^AT/,],  (55) 

n=l  n>m 

where  Ax/,  is  the  inner  hydrophone  time  delay  required  to  steer  a  beam  in  direction  0/,,  measured 
with  respect  to  the  receiver  array  axis. 


7.5.2  Flow  Noise 

Flow  noise  (Bradley  and  Wagstaff  1991)  in  RANDI  3.1  is  based  on  empirical  data.  Flow  noise 
is  assumed  to  produce  a  cross-spectral  density  matrix  Cpiow,mn  for  the  receiver  array  of  the  form 


C Flow, mn  —  L Flow  (,flfd)^mn  j 


(56) 
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where /is  the  frequency,/^  is  the  design  frequency  of  the  receiver  array,  Lpiow  is  the  omnidirectional 
flow  noise  at  the  hydrophone  level,  and  bmn  is  the  Kronecker  delta  function. 

Beam  output  levels  for  either  the  FFT  or  DFT  beamformers  due  to  flow  noise  are  independent 
of  steering  direction  and  are  specified  by 


N 

^  Flow,  b  ~  LpioMlfd)  2  \  •  (57) 

The  hydrophone  flow  noise  level  Lpiow  depends  upon  the  frequency  and  the  tow  speed.  It  is 
interpolated  from  the  empirical  data  for  frequencies  in  the  range  of  10  to  340  Hz  and  tow  speeds 
in  the  range  5  to  18  kt.  The  flow  noise  formulation  is  not  valid  for  Nflfd<  1. 

At  a  frequency  of  30  Hz,  an  array  with  64  elements  spaced  15  m  apart  results  in  omnidirectional 
flow  noise  levels  as  shown  in  Fig.  19  for  the  allowable  range  of  tow  speeds.  At  high  tow  speeds,  flow 
noise  can  easily  overwhelm  the  receiver  array. 


7.6  Combined  Array  Response 

Shipping,  wind,  system,  and  flow  noise  are  effectively  independent  sources  of  noise.  By  incoherently 
combining  model  predictions  of  the  array  response  to  these  four  noise  sources,  the  total  response 
can  be  estimated.  For  both  DFT  and  FFT  beamformers,  this  combination  is  performed  at  the  beam 
output  level.  The  result  for  a  particular  beam  is  simply 

^Arrayyb  =  Bship,b  +  ^Windyb  +  ^System, b  +  ^Flow^b  • 

The  incoherent  DFT  beamformer  was  used  to  bearaform  the  received  noise  from  the  East  China 
Sea  example  described  earlier.  Figures  20,  21,  and  22  give  the  beam  noise  results  from  the  rigorous, 
radial  accumulation,  and  high-resolution  radial  accumulation  shipping  propagation  methods, 
respectively.  The  three  results  are  plotted  on  the  same  scale  for  ease  of  comparison.  The  beam  noise 
includes  the  shipping  noise  field  (shown  in  Figs.  9,  10,  and  11),  as  well  as  noise  from  a  10-kt  wind. 
Beam  noise  results  using  the  radial  methods  compare  very  well  with  results  from  the  rigorous 
method.  The  high-resolution  radial  accumulation  method  is  only  slightly  different  from  the  radial 
accumulation  method. 


Fig.  19  —  Example  of  omnidirectional  flow  noise  levels 
vs.  tow-ship  speed  for  30  Hz  using  a  64-element 
horizontal  array  with  15-m  spacing 
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Fig.  20  —  DFT  beam  noise  with  incoherent  processing 
for  the  East  China  Sea  example  using  the  rigorous 
propagation  method  for  ships 


Fig.  21  — DFT  beam  noise  with  incoherent  processing 
for  the  East  China  Sea  example  using  the  radial 
accumulation  propagation  method  for  ships 
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Fig.  22  —  DFT  beam  noise  with  incoherent  processing 
for  the  East  China  Sea  example  using  the  high-resolution 
radial  accumulation  propagation  method  for  ships 


Fig.  23  —  DFT  beam  noise  with  coherent  processing 
for  the  East  China  Sea  example  using  the  high-resolution 
radial  accumulation  propagation  method  for  ships 


The  coherent  beam  noise  for  the  rigorous  propagation  method  is  shown  in  Fig.  23.  Beam  noise 
levels  near  endfire  are  somewhat  suppressed  when  coherent  processing  is  used  instead  of  incoherent 
processing  (compare  to  Fig.  20).  Coherent  processing  may  not  be  appropriate  for  the  radial  propagation 
methods  since  phase  information  is  either  not  included  or  only  approximated  for  individual  ships. 

The  coherent  and  incoherent  omnidirectional  levels  are  calculated  for  each  hydrophone  in  the 
array.  The  omnidirectional  level  for  the  array  as  a  whole  is  taken  to  be  the  median  of  the  hydro¬ 
phone  levels  in  decibels.  The  incoherent  omnidirectional  level  compares  favorably  for  the  different 
propagation  methods:  89.7  dB  for  the  rigorous  method  versus  88.7  dB  for  the  radial  methods.  The 
coherent  omni-directional  levels  are  equivalent  for  the  rigorous  and  radial  accumulation  methods, 
and  are  equal  to  89.1  dB. 
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